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CsJ ■ Abstract. We have studied the mass distribution in the lensing galaxy 

2237+0305 using constraints from both gravitational lensing and pilo- 
se ■ tometric and spectroscopic observations. We find that with sufficient 

dynamical information we can constrain the mass-to-light ratio of the 
Q^ ■ luminous components and determine the dynamical contribution of the 

disk. In addition, future observations should allow us to place constraints 
on the shape of the dark matter halo and its inner slope. 
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. 1. Introduction 

o . 

Gravitational lensing has become an important tool in recent years for studying 
many aspects of the universe. As it probes total mass, independent of the light 
distribution, it opens new windows onto studies of large scale structures where 
dark matter has dynamical importance. 

With new galaxy lenses being discovered each year, the statistical use of 
systems to study dark matter halo shapes and orientations is fast becoming 
possible. Most lensing systems lie at high redshift as this is where the lens- 
ing cross-section is optimal, however for detailed studies of the lensing galaxies 
themselves, lower redshift systems are better laboratories. 

2237+0305 was discovered as a low redshift quadruple lens by Huchra et 
al. (1985). At z=0.039 it is a visibly extended early-type spiral galaxy with 
multiple luminous stellar components. In addition, its proximity to us places 
its four lensed quasar images close to the centre of the galaxy at ~700pc (most 
lens images sit ~5kpc from the deflector's centre) making the probed region of 
the galaxy more compact. These features of the 2237+0305 system make it an 
ideal laboratory for studying the inner structure of the dark matter halo and 
properties of the lensing galaxy. 

Most galaxy lenses are early-type and studies of these systems can also 
provide important constraints on the dark matter (for example Treu, this meet- 
ing). Studies with spiral lenses are less frequent due to the fewer number known 
(elliptical galaxies present a larger lensing cross-section probability) but can pro- 
vide information about the luminous components of the galaxy (for example the 
system B 1600+434, Mailer et al. 2000). 

The importance of the stellar disk to the dynamical support of the galaxy 
is a controversial question in modern galactic astronomy. Sackett (1997) defined 
a maximal disk as one which provides 75-95 per cent of the galaxy's dynami- 
cal support. Whether a disk is maximal or sub- maximal has an impact on our 
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understanding of disk dynamics and is therefore an important quantity to deter- 
mine. Similarly, the mass-to-light ratios in the stellar components of a late-type 
galaxy are not well determined and their values can place constraints on star 
formation models. Traditional studies of spiral galaxies have used light profiles 
of the luminous components and spectroscopic dynamical information to deter- 
mine these properties of the galaxy. In a complex system with multiple mass 
components, a degeneracy exists between the mass-to-light ratio of the major 
stellar component and the mass scaling of the dark matter halo. The additional 
information provided by lensing can break this degeneracy. 

2. Galaxy Modelling 

2237+0305 has been studied extensively in the twenty years since its discovery 
(for example Yee 1988, Foltz et al. 1992, Schmidt 1996). Many simplified 
models of the galaxy's mass distribution have been published which focus on 
different aspects of the system. We aim to model the major mass components 
with physically reasonable profiles and use these models to constrain the mass- 
to-light ratios, disk maximality and dark matter halo inner slope. 

Schmidt (1996) was the first to investigate the effect of the galaxy's bar on 
the lensed image positions. He found the bar provides a significant torque to 
their positions and is thus a critical component to any lens model. In addition to 
the bar, we model an inclined stellar disk, an elliptical bulge at the same position 
angle, and a dark matter halo. The disk is modelled as a thin structure with an 
exponential surface density profile, the bulge with a de Vaucouleurs profile and 
the bar as a Ferrers ellipse, in line with previous works. The dark matter halo we 
model with both a flattened NFW profile (Navarro, Frenk &z White 1996) and 
as a softened isothermal sphere (SIS). The parameters characteristic to these 
models have been taken from previous observational studies and are presented 
in Table 1. The source position remains unknown and increases the number of 
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Table 1. Fixed and free parameters and their source for the mass 
models of the four major components. The label 'free' refers to free 
parameters in the model and r c and denote the core radius and scale 
length respectively. 1 Schmidt 1996, 2 Yee 1988, 3 Racine 1991, 4 Irwin et 
al. 1989. 

free parameters by two. In addition to the constraints from observations, the 
four image positions are well determined by HST (positions of Blanton, Turner 
& Wambsganss 1998 used in this study). Barnes et al. (1999) observed the 



APS Conf. Ser. Style 



3 



neutral hydrogen 21cm emission line from the galaxy with the VLA and this 
low resolution data has allowed us to obtain two independent points in the outer 
rotation curve. 

With the known image positions and rotation curve information, there are 
ten constraints on the model. For the softened isothermal sphere, there are 
seven free parameters and for the flattened NFW there are eight. The problem 
is therefore over-constrained. 

3. Inversion Method 

The unknown parameters are varied to fit both the image positions and rota- 
tion curve points using a basic minimisation procedure and the \ 2 parameter. 
The disk and bulge scale lengths are also allowed to move within their la un- 
certainties in order to explore the parameter space adequately. The bending 
angles for each mass component are computed and their relative contributions 
(mass-to-light ratios) varied to fit the constraints. Solutions with unphysically 
large or small values for the mass-to-light ratios were discounted. Since most 
components contribute shear as well as convergence to the lensing configuration, 
degeneracies exist between some parameters. 

4. Results 

Statistically acceptable solutions were found for both dark matter halo density 
profiles. The degree of ellipticity required for the flattened NFW profile was very 
small and consistent with zero. Both solutions provided similar convergence 
within the image regions and therefore allowed the same combination of the 
luminous mass components. 

For 3 (2) degrees of freedom, an acceptable x 2 value at la confidence is 3.6 
(2.4). Table 2 shows the best-fitting parameters for the models. It is evident 
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0.01 


Scale 


r c = 1.03±0.02kpc 


r h = 12.0±0.6kpc 


x 2 


1.4 
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Table 2. Best fitting parameter values for the two halo density pro- 
files. Note the similarity between mass-to-light ratio values and the 
almost spherical NFW halo. 

from the similarity of the mass-to- light ratios and the small NFW flattening that 
these two halo profiles would contribute comparable convergence and shear to 
the lensing. Figure 1 shows the rotation curve of a galaxy with the best-fitting 
parameters for the softened isothermal sphere halo and luminous components. 
The bulge component clearly dominates the centre of the galaxy, providing the 
majority of the dynamical support within the image region. In the outer regions 
of the galaxy the dark matter component dominates, as expected from results of 
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Figure 1. Best-fitting rotation curve based upon rotation points from 
HI observations and image positions for a softened isothermal sphere 
dark matter halo. The solid line is the total rotation curve, the dash- 
dotted line is the dark matter halo, the dashed line is the bulge, the 
dotted is the disk and the long-dashed is the bar. 

previous work deconvolving galactic rotation curves. The best-fitting rotation 
curve using the flattened NFW dark matter halo is essentially the same. 

The disk is clearly sub-maximal according to the definition of Sackett (1997), 
however the importance of the bulge in this system should not be forgotten. This 
early-type spiral is not the ideal system for studying the importance of disk stel- 
lar support as the bulge clearly dominates the disk. 

5. Discussion 

Although the results presented are acceptable fits to the data, they are not likely 
to describe the real situation in the galaxy. Firstly, the flux ratios of the images 
have not been included in the minimisation procedure, even though they are 
additional constraints. The first reason for this is the uncertain nature of ra- 
tios measured for most lensing systems. Q2237+0305 is not a radio-loud quasar 
and radio measurements are difficult. Optical fluxes probe a small region in 
the source plane and are subject to possible microlensing leading to spurious 
values. Agol, Jones & Blaes (2000) measured mid-IR fluxes for this system ar- 
guing the larger extent of the emitting region would reduce the contamination 
of microlensing and the waveband would minimise the effects of dust redden- 
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ing. We will use these measurements to compare with our best-fitting model. 
The second reason for omitting the flux ratios in the analysis is the contention 
over their applicability when dark matter substructure is taken into account 
(Schechter & Wambsganss, this meeting). Our dark matter is modelled as a 
smooth component. 

Comparison of the flux ratios shows excess flux in one of the images with 
both profiles for the dark matter. This result is quite robust and suggests an 
underlying problem in the mass distribution. 

In addition to this problem, the bulge mass-to-light ratios are extremely low 
compared with other studies of the bulge component in I-band (see for example 
Fukugita, Hogan & Peebles 1998) which place the figure closer to M/L=5.0T. 
Clearly, increasing the bulge mass changes the convergence and the shear in the 
image region due to its dominance of the mass there. Without changing the 
remainder of the mass components, the image positions would be torqued out 
of place. 

The location of the images defines the projected mass within that region 
quite accurately (Wambsganss & Paczynski 1994 measured M enc = (1.48±0.01) 
x lO lo M h7 5 1 ) and therefore the bulge mass-todight ratio can only be increased 
so far before it alone provides all of the convergence. In addition, as the bulge 
is an elliptical component, increasing its contribution will increase the shear. 

The one completely unknown component of the mass is the dark matter 
halo. The best-fitting disk and bar sizes are what one would expect. If the bulge 
is too small in our models then it would appear the dark matter distribution 
would need to be changed to accommodate this. The most likely outcome is a 
halo with little contribution within the image region but sufficiently flattened as 
to counteract the increased torque from the bulge. Without further information, 
this is purely speculation. 

Additional dynamical information can help to answer these remaining ques- 
tions. With only rotation data for the outer galaxy, critial regions of the rotation 
curve are unknown. Well measured points in the galactic centre will tightly con- 
strain the allowed dynamics in the important image region, as shown in Figure 
2. It compares the expected observed rotation curve for this galaxy with a max- 
imal bulge and the two different halo profiles. This solution does not fit the 
image positions, but shows the possible difference in the inner rotation curve. It 
is this splitting of curves that we hope to probe with future spectroscopy of the 
galaxy. 

6. Conclusions 

We have modelled the galaxy 2237+0305 with multiple mass components to 
investigate the mass-to-light ratios of the luminous components, the dynamical 
importance of the disk and bulge and the shape of the dark matter halo. We 
require additional dynamical information to define the dark matter contribution 
accurately. This system is potentially useful for constraining the mass-to-light 
ratio of the bulge component and the shape and inner slope of the dark matter 
halo. 
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